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Abstract 

A wide class of Higgs sectors is investigated in supersymmetric standard models. When the 
lightest Higgs boson (h) looks the standard model one, the mass (rrih) and the triple Higgs boson 
coupling (the hhh coupling) are evaluated at the one-loop level in each model. While is at 
most 120-130 GeV in the minimal supersymmetric standard model (MSSM), that in models with 
an additional neutral singlet or triplet fields can be much larger. The hhh coupling can also be 
sensitive to the models: while in the MSSM the deviation from the standard model prediction is not 
significant, that can be 30-60 % in some models such as the MSSM with the additional singlet or 
with extra doublets and charged singlets. These models are motivated by specific physics problems 
like the //-problem, the neutrino mass, the scalar dark matter and so on. Therefore, when h is 
found at the CERN Large Hadron Collider, we can classify supersymmetric models by measuring 
rrih and the hhh coupling accurately at future collider experiments. 

PACS numbers: 14.80.Da, 12.60.-i, 12.60.Fr 
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I. INTRODUCTION 



Physics of electroweak symmetry breaking (EWSB) is the last unknown part in the stan- 
dard model (SM). Experimental identification of the Higgs sector has been one of the most 
important issues in high energy physics. Yet no Higgs boson has been found at the CERN 
LEP and the Fermilab Tevatron, the data are used to constrain the Higgs boson parame- 
ters jl, 2]. The Higgs boson search has also started at the CERN Large Hadron Collider 
(LHC). When a scalar boson is found in near future, its property will be measured as ac- 
curately as possible to see whether the scalar is really the Higgs boson. In particular, to 
understand the essence of EWSB, we must know the self-coupling in addition to the Higgs 
boson mass. The triple Higgs boson coupling may be explored at the LHC and its upgraded 
version^, the International Linear Collider (ILC)|4] and its 77 option[5] or the Compact 
Linear Collider (CLIC). 

Experimental survey for the Higgs sector is important not only to confirm our picture 
for EWSB, but also to explore new physics beyond the SM (BSM). In the SM, appearance 
of the quadratic ultraviolet divergence in the one-loop calculation of the Higgs boson mass 
is a serious problem, which is called the hierarchy problem. This has to be eliminated in 
a model of BSM. On the other hand, we already know several phenomena of BSM such as 
dark matter, tiny neutrino masses, and baryon asymmetry of the universe. 

Supersymmetry (SUSY) is an attractive idea as a new physics scenario at the TeV scale. 
First of all, the hierarchy problem can be solved: i.e., the quadratic divergence due to 
a particle in the loop is cancelled by that due to the super partner particles. Second, a 
supersymmetric standard model (SSM) can naturally contain the candidate for dark matter; 
i.e., the lightest SUSY particle, whose stability may be guaranteed by the R-parity. One of 
the important predictions in the minimal SSM (MSSM), where the Higgs sector is composed 
of two Higgs doublets, is to the mass (ra/J of the lightest CP-even Higgs boson (h), which 
is predicted to be less than the Z boson mass (m z ) at the tree level, but can be maximally 
30-40 GeV above mz by radiative effects of the top quark and its super partners (stops) [6(. 
This prediction meets the present LEP bound on mn If such a light h will not be found, 
the MSSM must be ruled out. Is SUSY itself ruled out then? The answer is "No". In 
fact, there can be several variations for the Higgs sector in SSMs. A simple extension may 
be addition of a neutral gauge singlet field to the MSSM, which is known as the Next-to 
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MSSM (NMSSM) jz-jj. Solving the fi problem 10| may be a motivation to the NMSSM. 
It is known that in the NMSSM can be higher than that in the MSSM because of the 
additional contribution to rrih from the tri-linear term of the Higgs doublets and the singlet 
in the superpotential. In addition to the MSSM and the NMSSM, a model with additional 
scalar boson with higher representation can also push up rrih 11] ■ Furthermore, models with 
extended gauge sector can also have higher m h than the MSSM prediction |12|, depending 
on the model parameters. Therefore, rrih can be an important tool to discriminate SSMs. 

In this Letter, we investigate Higgs sectors in a wide class of SSM, when only the lightest 
Higgs boson h appears at the electroweak (EW) scale and the other additional particles such 
as extra Higgs bosons and SUSY partner particles are rather heavier (but not too heavy). 
The coupling constants of h to the SM particles are then similar to those in the SM at the 
tree level. In each model, rrih as well as the deviation in the triple Higgs boson coupling (the 
hhh coupling) are evaluated at the one-loop level. Their possible allowed values as well as 
their correlation are studied under the theoretical condition and the current experimental 
data. We here consider not only the MSSM and the NMSSM, but also further possible 
extensions of the MSSM. In order to widely examine various possibilities, we here do not 
necessarily require the paradigm of the grand unification. Instead, we may even retain the 
possibility of appearance of strong dynamics at the TeV scale as discussed in Ref. [3]. 



II. EXTENDED SUSY MODELS 



One way of the extension of the MSSM may be adding new chiral superfields such as 
isospin singlets (neutral S, singly charged Q± or doubly charged K±), doublets (H' u and 
H' d ), or triplets (£ with the hypercharge Y = or \± with Y = ±1), whose properties are 
defined in Table HI As we are interested in the variation in the Higgs sector, these new 
fields are supposed to be colour singlet. Although there can be further possibilities such as 
introduction of new vector superfields which contain gauge fields for extra gauge symmetries, 
models with extra dimensions, those with the R-parity violation, etc, we here do not discuss 
them. Eventually, in addition to the MSSM, we consider nine models in Table [Til In the 
analysis of this letter, triple coupling terms with doublet Higgs superfields and additional 
chiral superfields play an important role. In Table [TTJ such relevant terms are also listed. 
For anomaly cancellation, charged superfields are introduced in pair in each model. 
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TABLE I: Properties of additional chiral superfields. 
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W D X -^H U ■ X -H u + X -^H d ■ X+ H d 


Model-6 




• • 






Model-7 


• 


• • 




W D X H u H d S H u ■ H d S + X H' u H d S H u ■ H d S 
+ X H u H' d S H u ' H 'd S + X H' u H' d S H 'u ' H 'd S 


Model-8 




• • 
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W D X H u H d ( H u ■ Z H d + X H' u H d £_ H 'u ■ £ H d 
+ X HM H u ■ i H 'd + X H'M H 'u ■ ^ H 'd 


Model-9 


• • 


• • 




W D X HHn- H u ■ H'J 1 - + X HHU+ H d • H d^+ 



TABLE II: Particle entries in each SSM. Additional terms in the superpotential which are relevant 
to the extra chiral superfields are also shown in each model. 



These extensions can be motivated by solving various physics problems. For example, 
models with additional charged singlet fields can be used for radiative neutrino mass genera- 

l_l O O O 

tion [14j. Those with additional doublet fields may be required for dark doublet models 15], 



and the model with triplets may be motivated for the SUSY left-right model[16| or those 
with so-called the type-II seesaw mechanism [ttI ]. Although models in Table ITT1 can be im- 
posed additional exact or softly-broken discrete symmetries for various reasons, we here do 
not specify them as they do not affect our discussions. 
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III. THE METHOD 



The effective potential for the order parameter ip can be written at the one-loop order 
as(l8| 



(-l) 2s fN Cf N Sf 



f 



64vr 2 



In 



Q 2 



where /i 2 , and A are the bare squared mass and the coupling constant, m,f(ip) is the field 
dependent mass of the field / in the loop, and Nq, and Ng. are the degree of freedom of 
the colour and the spin with Sf being the spin of the field /. In SSMs, the Higgs sector has 
a multi-Higgs structure. When the extra Higgs scalars are heavy enough, only the lightest 
Higgs boson h stays at the EW scale, and behaves as the SM-like one. The effective potential 
in Eq. ([TJ) can then be applied with a good approximation. The vacuum, the mass m^ and 
the hhh coupling constant Xhhh are determined at the one-loop order by the conditions; 
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(2) 



where v (~ 246 GeV) is the vacuum expectation value for the EWSB. 



IV. THE MASS AND THE TRIPLE COUPLING OF THE LIGHTEST HIGGS 
BOSONS 

A. The Mass of the lightest Higgs boson 

In the MSSM, rrih is evaluated by using Eq. (J2]) as 

m\ = m| cos 2 2(3 + Smf oop , (3) 
where tan/3 = (H®)/ (H®). The first term of the R.H.S. comes from the D-term at the tree 



Q to satisfy the 



level. The quantum correction Sm 2 has turned out to be important {(3 
LEP data (m h > 114 GeV) |l|; 

8m 2 ~ ^ In ?M + 3m * X ? Sin2 ^ In < + OfX'/m* ) (4) 

where m 4 is the top quark mass, are the masses of stops f x and t 2 < m t 2 )> an d -^Q 
is defined such that the coupling constant of ti-t^-h is given by X 4 sin (3 / Because the 
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TABLE III: Formulae for m h and Xf^/Xlfh in various SSMs at the one-loop order [2jJ. Masses 
of the component fields S c , Tf, and flj are given by ms c , tut? and rasi,, while M| c = m 2 c — 
\ 2 HH( f ) .v 2 Ci/2, where 0£ = S c , T- c , or Oj and C_ = sin 2 /? and C+ = cos 2 /3 in both Model-5 and 
Model-9 and Ci = 1 in the other models. In Model-2, Model-7 and Model-8, there are similar 
non-decoupling effects to Model- 1. 



top Yukawa coupling is of order one, Smf oop can push up the upper bound to 120-130GeV 
when stop masses are of order one TeV. Higher order calculations for rrih have been studied 
in literature [201] . 

In general cases with additional chiral superfields, the F-term can also contribute to m^ 
at the tree level. For example, in Model- 1 (i.e., the NMSSM), the term of \hhsH u ■ HdS in 
the superpotential gives the coupling of X 2 HHS (H U ■ Ha) 2 in the Higgs potential, which yields 



m 2 h ~ m 2 z cos 2 2/3 + (X 2 HHS v 2 /2) sin 2 2/3 + 5m 2 oop (Model-l), 



(5) 



where 



^ m ioop( M odel-l) = <5m loop + 
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where m Sl , m Sa , and m- are masses of a CP-even scalar component, a CP-odd scalar 




FIG. 1: The upper bounds on rrih as a function of tan/3 for fixed values of Xrhs and X HHx _ = 
^HHx+ = ^ HH x 111 Model- 1 and Model-5, respectively. The red-filled region indicates the possible 
allowed region in the MSSM. 

component and a fermion component of S, respectively. We here assume that the mixing 
between c omp onent fields are so small that each component field can be regarded as a mass 
eigenstate 



301 ] . We consider the case that the scalar and the fermion components of the 
additional fields are approximately degenerate and the logarithmic correction from them are 
negligible. Due to the X 2 hhs t errr b m h can be above 130 GeV. It is bounded from above by 
the renormalization group equation analysis assuming the condition of avoiding the Landau 



pole below a given value of A [8|, |22| . The upper bound of mh is evaluated to be about 140 
GeV at tan/3 ~ 2 (about 450 GeV at tan (3 ~ 1) for A ~ 10 16 GeV (4 TeV) ^. 

Similarly, the enhancement of due to the F-term contribution can also be realized in 
models with an additional field <fi where the superpotential has the gauge singlet operator 
like H u H d (j>. In Table HQ Model-1, Model-2, Model-5, Model-7 and Model-8 satisfy this 
condition, where mh receives the F-term contribution and can be significantly larger than 
that in the MSSM. Approximate formulae for m/, are given in Table IIHI 

In Fig. [H the upper bounds on in Model-1 and Model-5 are shown as a function of 



tan ft, and the possible allowed region in the MSSM is also indicated by the red-filled region. 
The coupling constants Xhh4> (4> — S and x±) are taken as < Xhh4> < 2.5, whose upper 
limit corresponds to A > 4 TeV in Model- 1 for the case with the Topical scale of the soft 
SUSY breaking to be 500 GeV[3l| . The detail is shown elsewhere 21J. In Model-1 with a 
fixed value of X HHS , mu can be maximal for tan ft = 1, while in Model-5 it becomes maximal 
for large values of tan. ft for a fixed value of Xhh\- The maximal value in Model-1 becomes 
asymptotically the same as that in the MSSM in the large tan/3 limit up to the one-loop 
logarithmic contributions. 



B. The hhh coupling 



We turn to the discussion on the quantum effect on the hhh coupling in the SSMs in 
the case where h is regarded as the SM like Higgs boson. To this end, we start from the 
case in the non-SUSY extended Higgs sector. It is known that in the non-SUSY two Higgs 
doublet model (THDM), the hhh coupling can receive large non-decoupling effects from 
the loop contribution of extra Higgs bosons, when their masses are generated mainly by 



EWSB 



23 



24j. When h is the SM like Higgs boson, physical masses of the extra scalar 
bosons are expressed by 



M 2 + 



kv 2 



(7) 



where <3>j represents H°, H ± or A , and M is the invariant mass scale which is irrelevant to 
the EWSB, and A is a coupling for $>\<& i hh. The physical meaning of M is discussed in, 



for example, Ref. 



calculated as 
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24j . The one- loop contribution to the renormalized hhh coupling is 
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One finds that for M 2 3> X^v 2 it becomes 

,2 
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(9) 



which vanishes in the large M limit according to the Appelquist-Carazzone decoupling 
theorem {25 1 . On the contrary, when the physical scalar masses are mainly determined by 
the X { v 2 term, the loop contribution to the hhh coupling does not decouple, and the quartic 
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FIG. 2: Possible allowed regions in the m/ 1 -(AA/ l / l / l /A/ l / l / l ) plane in the MSSM, Model-1, Model-5 
and Model-9 for each tan/3 value. We scan the parameter space as < Xhh^ < 2.5, 0.5 TeV < 



TO r < 1.5 TeV, and 0.5 TeV < rrii C for each model. 

£1,2 fi 



powerlike contributions of m $ remain; 



\THDM 

A hhh 
\SM 
A hhh 



1 + 



(m%o + m\ + 2m A H ±) 



(10) 



Consequently, a significant quantum effect can be realized for the hhh coupling when m| a > 



in 



;. The size of the correction from the SM value can be of 100% for 



120GeV, M ~ 0, 



and m H „ ~ 



rriAo 



rriH± — 400GeV under the constraint from perturbative unitarity 



26). 



Such a large non-decoupling effect on the hhh coupling is known to be related to the strongly 



first order EW phase transition [27J which is required for the EW baryogenesis 



28| 



Let us discuss the hhh coupling in the SSMs listed in Table |TTJ In the MSSM, it is 
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evaluated at the one-loop level as 

4 ( rriiiml + m 2 ) 3 Y 2 ,, 2 cin 2 R } 

(11) 
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up to the stop mixing 271 ] . This result coincides with that in the SM in the decoupling limit 



24j. On the contrary, in a class of the SSMs; Model- 1, Model-5 and Model- 



9, the non-decoupling effect can appear in the hhh coupling, similarly to the case of the 
THDM. Approximate formulae of the A^ del /Af^ are given in Table IIII1 In the following, 
we explain the results in each model in order. 

In Model- 1 a scalar boson Si and a pseudo-scalar boson S2 from extra isospin singlet 
field S are running in the one-loop diagrams of the hhh coupling. Their physical masses 
m| (c = 1, 2) are given by 



< = M 2 Sc + -\ 2 HHS v 2 , (12) 



1 

2' 

where M$ c represent the invariant mass parameters and Xhhs is a coupling constant of the 



H u ■ H^S in the superpotential. The loop effect decouples in the large M| limit in the same 
way as the THDM. The non-decoupling property appears when \ 2 HHS v 2 > M| . However, 
differently from the THDM, M\ may not be too small, because this is directly related to the 
mass of the singlino. When Mg c is taken to be around 500 GeV, the one-loop contribution 
still turns out to be important. For example if M$ c — 500 GeV, Xhhs — 2.5 and tan/3 ~ 20 
(ms c — 660 GeV), the correction to the SM prediction can be as large as 40%. Here we take 
nig ~ ms- The quantum correction to the hhh coupling in Model- 1 strongly depends on 
tan/3 because the lightest Higgs boson mass rrih depends on tan/3 (see Fig. [3]). 

As is seen in Table 1111} similar effect can also be realized in Model-5 and Model-9. In 
Model-5 there are two types of the triplet fields (T_ and T + ), each of which gives six degrees 
of freedom, namely two neutral scalar bosons, a pair of singly charged Higgs bosons, and a 
pair of doubly charged Higgs bosons. However, only neutral and singly charged degrees of 
freedom contribute to the loop effect of the hhh coupling, because a term like hhT ++ T 
cannot exist in the Higgs potential. In this formula, for each i (= — or +), T} is the CP- 
even scalar boson, Tf is the pseudo-scalar boson, and (Tf ± iTf) / y/2 represents a pair of 
the singly charged Higgs bosons. Their physical masses m^c are 

m 2 Tr = M 2 . + ^X HHxi v 2 Ci , (13) 
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where Myp are the invariant mass parameter, and C_ = sin 2 [3 and C + = cos 2 /3. Here 
m T 3 = m T A (= m T ±) and M T 3 = M t a are assumed. One might think that similarly to 
the THDM and Model-1, the correction to the hhh coupling can be significant, when mjp 
is mainly from the \hh Xi v 2 term. However, this is not the case because rrih becomes also 
large by the same mechanism of enhancement as the hhh coupling so that the net correction 
cannot be very significant. 

Next we discuss the case of Model-9 where there are two extra doublets and two singly 
charged singlets. Notice that the extra doublets do not contribute to the one loop correction 
to the hhh coupling because there is no corresponding F-term, so that only the one-loop 
effect of charged singlets can be important. Two singly charged scalar bosons with i = — 
or +) are running in the one-loop diagram of the hhh coupling. Their physical masses tuq. 
are given by 

ml, = K t + l^nn^a , (14) 

where are the invariant mass parameters. Differently from Model-5, we can expect large 
deviation in the hhh coupling from the SM prediction in this model, because rrih does not 
get a significant enhancement from the F-term contribution. We have similar large effect to 
that in Model-1 and THDM, when mo is mainly from AlUro v 2 - 

* i 

In summary, a non-decoupling quartic power-like scalar mass effect occurs in the hhh 
coupling when additional scalar field (such as Si )2 , H' d , 7±' + , etc.) receives its mass mainly 
from v via the operators HiHj(j)(j), which are generated from HiHjip in the superpotential, 
where Hi (i = u, d) are the Higgs doublets. This effect appears in Model-1, Model-2, Model- 
5, Model-7, Model-8 and Model-9, while there is no such contribution in the other models. 
Even in models with this effect, the deviation from the SM prediction in the hhh coupling 
cannot be significant, when rrih is also enhanced by the same F-term contribution. For 
example, in Model-1 with a small tan/3 (~ 1), large X HHS gives a significant contribution 
to the numerator. However at the same time, the m h in the denominator gets the large 
contribution as seen in Eq. Consequently, enhancement by the non-decoupling effect 
of particle S are weakened, and the deviation in the hhh coupling is not very important. 
Model-2, Model-5, Model-7 and Model-8 also correspond to this case. On the other hand, in 
Model-9 for example, the F-term of the operator H U -H' U Q^ or Hd-H' d Q + cannot contribute to 
rrih but gives a quartic power contribution of rriQ +/ _ in the one-loop corrected hhh coupling, 
so that the deviation in the hhh coupling from the SM value can be significant. 
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C. The correlation between m h and the hhh coupling 

We scan the parameter space in each model to find allowed regions in the 77^- 
(AAj^f el /A™ h ) plane under the assumption of \hh^> < 2.5 at the EW scale, where 
^■^hhh Cl = ^hhh d ~ ^Thh- I* 1 Fig- El we show the possible allowed region for several value 
of tan/3 = 1, 4, 8 and 20. The coupling constants Xhhq (0 = S, x± an d ^±) are taken to 
be less than 2.5 as in Fig. [TJ The stop masses are scanned as 0.5 TeV < tti^ < 1.5 TeV. 
We also scan the physical masses of the extra scalar bosons as 0.5 TeV < m,. The mass of 
fermion component is taken as same as the mass of the scalar component for each extra field. 
We note that the parameters are scanned such that the additional contributions to the rho 



parameter are negligible [32J. The region in the MSSM is indicated as the red-filled one. The 
possible allowed region in Model- 1 depends largely on tan j3: for smaller (larger) tan /3, 
can be higher (lower) and AAj^ el_1 / is smaller (larger). Model-5 is relatively insensitive 
to the value of tan/3: can always be larger than about 300 GeV while AA^ del_5 /A^ 
remains less than about 10 %. On the other hand, in Model-9, although the possible value of 
rrih is similar to that in the MSSM, the deviation in the hhh coupling can be very large: i.e., 
^^hhh Cl ~ 9 '/ '^hhh ~ 30 — 60 %[33]. When we consider the higher value of A, which corresponds 
to the smaller upper bound on X HH ^, the possible allowed region becomes the smaller. 

In Fig. [31 possible allowed regions with scanned tan (3 are shown in the m/ l -(AA^f e1 / A|jjJ 
plane in Model- 1, Model-5 and Model-9 as well as the MSSM. The maximal values of Xjjjj^ 
in Model- 1, Model-5 and Model-9 are taken to be the same as those in Fig. |2j The region 

is 



in the MSSM (Model- 1 with < X H hs < 0.75, which corresponds to A ~ 10 16 GeV 
indicated as the red-filled (cyan-filled) one. The possible allowed regions are different among 

Mode 
hhh 



the models so that the information of and AAj^f 61 can be used to classify the SSMs 



V. DISCUSSION AND CONCLUSION 



We have studied mn and the hhh coupling at the one-loop level in various SSMs, where h 
is the lightest SM like Higgs boson. In a class of SSMs, the mass of the lightest Higgs boson, 
which is lower than 120-130 GeV in the MSSM, can be much higher due to the F-term 
contribution. Such an enhancement appears in the models with the extra neutral singlet or 
the triplet superfield. The upper bounds on nih are determined by the size of the coupling 
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FIG. 3: Possible allowed regions in the rrih-(A\hhh/ Xhhh) plane in the MSSM, Model-1, Model-5 
and Model-9 with scanned tan f3. 

constants of the F-term, which can be constrained by the renormalization group equation 
analysis with an imposed cut-off scale A. Consequently, rrih can be higher than 300-400 GeV 
when A is at the TeV scale in Model-1 with small tan (3 values and Model-5. 

Although the one-loop correction to the hhh coupling due to the extra scalar components 
vanishes in the decoupling limit, it can be significant in particular SSMs such as Model-1 with 
large tan/3 and Model-9, when A^^.f 2 ~ M|. where is the invariant mass parameter of 
the extra field 0j in the loop. In such a case, quartic powerlike mass contributions can appear 
as non-decoupling effects, and the correction can be larger than several times ten percent 
under the constraint from parturbativity. In this letter the analysis has been restricted 
in the effective potential method, where all the external momenta are set on zero. In the 
actual measurement of the hhh coupling, one might think that the momentum dependences 
would be important. For example, at the LHC the hhh coupling may be measured by W 
fusion process yF +,| W~* — > h* — y /i/i|3j, where the measured hhh coupling is a function 
of the energy of the elementary process. At the ILC and its 77 option, the processes 
e + e" — > Z* — y Zh* — y Zhh[4] and 77 — y h* — y hh\^ can be used. The energy dependence 



of the hhh coupling has been discussed in Ref. [23J: see Fig. 3 in it. It is shown that unless 



\/I ^> 2m/,, the energy dependence is small in the bosonic loop contributions to the hhh 
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coupling |34j . In addition, the non-decoupling effect in the wave function renormalization is 
at most quadratic instead of quartic in mass. Hence the wave function correction to the hhh 
coupling is known to be as large as of order one percent and it is negligible (2^. Consequently 
the calculation by using the effective potential gives a good approximation for our analysis. 

In conclusion, even when only h is observed in future, precision measurements of rrih and 
the hhh coupling can help discriminate the SSMs. Such discrimination can be improved if 
extra information for tan can be used from, for example, future flavour experiments such 
as those for B s — > fifj,, B — > tv and so on. 

In any case, the hhh coupling is required to be measured with 0(10) % accuracy, which 



may be expected at future colliders such as the LHC upgrade, the ILC and the CLIC 



4]- 
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